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Abstract 

A recipe is given for creating material with a desired refraction coefficient by em- 
bedding many small particles in a given material. To implement practically this recipe, 
some technological problems are to be solved. These problems are formulated. 

PACS 43.20.+g, 62.40.+d, 78.20.-e. 
MSG 35J10, 74J25, 81U40 

Keywords: metamaterials; nanotechnology; refraction coefficient; negative refraction, 
wave focusing 

1 Introduction 

In a series of papers [6]-|13j and in the monograph [3] the author has developed wave 
scattering theory by many smah bodies of arbitrary shapes and, on this basis, proposed 
a method for creating materials with some desired properties. The goal of this paper is 
to explain this method and to formulate two nanotechnological engineering problems that 
should be solved in order that the proposed method can be implemented practically. 

Let us formulate the statement of the problem. Let D C be a bounded domain 
filled in by a material with a known refraction coefficient. The scalar wave scattering 
problem consists of finding the solution to the Helmholtz equation: 

Lqu := {V^ + k^nlix))u = in M^ Imno(x) > 0, (1) 



u = uq + v, Mo := e'''"'^', (2) 

V = Ao(/3, h o(-), r := Ixl — > oo, /?:=—. (3) 

The function no(x) is assumed bounded, piecewise - continuous, 

n{x) = 1 in D' := R^\D, lmnl{x) > 0. (4) 
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The function A{P,a) is called the scattering amplitude. The wavenumber k = where 
A is the wavelength in D', a G S'^ is the direction of the incident plane wave uq, is the 
unit sphere in R^, /3 G 5^ is the direction of the scattered wave. The solution to problem 
([T])-([3]) is called the scattering solution. It is well-known ([2j) that this solution exists and 
is unique under our assumptions, namely nQ(x) is bounded sup^jgi^s |ng(x)| < no = const, 
Imno(x) > 0, no(x) = 1 in D', so that AQ{f3,a) is determined uniquely if no(x) is given. 
We assume A: > fixed and do not show the dependence of AQ{P,a) on k. The operator 
Lq at a fixed A; > can be considered as a Schroedinger operator Lq = + A;^ — qo{x), 
where qo{x) := k'^nQ{x) — k"^. 

Problem(P): M^e want to construct a material in D with a desired refraction coeffi- 
cient, i.e., with a desired function n'^{x), Imn^(rE) > 0. 

Why is this problem of practical interest? We give just two reasons (out of, possibly, 
many). First, creating a refraction coefficient such that the corresponding material has 
negative refraction is of practical interest. One says that a material has negative refraction 
if the group velocity in this material is directed opposite to the phase velocity. 

Secondly, creating such a refraction coefficient that the corresponding scattering ampli- 
tude A[I3) = j4(/3, a) at a fixed a G 5^ (and a fixed A: > 0) approximates with any desired 
accuracy e > an arbitrary given function /(/?) € L?'{S'^) is of practical interest. This 
problem we call the problem of creating material with a desired wave-focusing property. 

In Section 2 we explain how to create such a material. In Section 3 we discuss the 
technological aspects of our recipe for creating material with a desired refraction coefficient. 
Creating material with negative refraction is discussed in Section 4. 

2 Creating material with a desired refraction coefficient 

The basic idea of our method is to embed in D many small particles Dm, 1 < m < M, with 
a prescribed boundary impedance. The number of the embedded particles in an arbitrary 
open subdomain A C D is given asymptotically as 



Here </> > if a > 0, 0(0) = 0, is strictly monotonically growing, N{x) > is a 
piecewise-continuous function which we can choose as we want. The smallness of the 
embedded particles means that A;a <C 1, where a = ^ max^ diam j where Dm is the 
m-th small particle. Let M = M{a) be the total number of the embedded particles. For 
simplicity we assume that Dm is a ball, centered at a point Xm G D, of radius a. The 
distance d := miiim^j dist{Dm, Dj) between neighboring particles is much larger than a, 
d ^ a, but d is not assumed much larger than A, so in our theory there can be many 
small particles on the wavelength in D. This means that the distribution of the embedded 
particles is not necessarily "diluted". 

The scattering problem in the region with embedded particles is formulated as follows: 




a 



0. 



(5) 




(6) 
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Un = CmU on Sm, l<m< M, (7) 
dV 1 

U(x) = u(x) + V, ikV = o(-), r^oo. (8) 

or r 

Here is the unit normal to Sm, pointing out of D^, Sm is the surface of D^, M is the 
total number of the embedded particles Dm, Cm is the boundary impedance, Im^m < 0, 
u is the scattering solution in the absence of the embedded particles, i.e., the solution to 
problem ([I])-©. 

Let us assume that 

Cm = — - — , — 1<J<<2, x = const, (9) 
where h{x) is a piecewise-continuous function in D, 

lm/i(x)<0. (10) 

The function N{x) > in (UD and h{x), Imh < 0, we can choose as we wish. These 
functions will define the desired refraction coefficient of the new material which is obtained 
in the limit a — > 0. We assume that 

M = M{a) = 0{a''-'^), d = 0{a^), (11) 
and the function (j){a) in ^ is 

(t){a) = a^-''. (12) 

The solution to problem ©-([T]) depends on a, C/ = U{x,a), but we do not show this 
dependence for brevity. We also do not show the dependence of C/ on q and /c, because 
they are fixed. Our standing assumptions are ([5]), ([9])- p^ . They are not repeated in 
the statements of Theorems. 

Denote by g{x, y) the Green's function of the operator Lq in the absence of the em- 
bedded particles, Log{x,y) = —6{x — y) in B? , g satisfies the radiation condition. 

Theorem 1 Problem ©-([Tl) has a solution, this solution is unique, and it is of the form: 

M „ 

U{x) = u{x)+^ / g{x,s)am{s)ds, (13) 

m=l 

where the function are uniquely determined by the data, i.e., by UQ{x,Oi), Cm; D^, 
1 <m < M , and n'^{x). 

Theorem 2 The solution U to problem ([6]) -dH) has the following asymptotic as a ^ 

M{a) 

U{x) = u{x) + g{x,Xm)Qm + o{l), a ^ 0, ini\x -Xm\>d, (14) 
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where 

Qm-= / cr„i{s)ds = -4-Kh{xm)ue{xm)a^^'^[^ + o{l)], a ^ 0, (15) 

J Sm 

CTmis) = -h{Xm)UeiXm)a~'^[l + o{l)], O ^ 0, (16) 

and 

Ue{x„,) :=u'i^\xm), u'f^^x) := u{x) + j g{x,s)am'{s)ds. (17) 

Theorem 3 Assume that f is an arbitrary bounded piecewise continuous function in D, 
and the points distributed in D so that ([5|) holds. Then 

\ixn 6{a) V /(x„) = / f{x)N{x)dx. (18) 

Theorem 4 Under the assumptions dH), (0), ([^- (|1'2|) the following limit exists: 

lim u^r\x) = Ue{x), Vx G D. (19) 

This limit does not depend on m and is the unique solution to the equation 

Ue{x)=u{x)+ / g{x,y)p{y)ue{y)dy, (20) 
Jd 

where 

p{x) = 47riV(x)/i(2;). (21) 

Theorem 5 The limiting material in D, which is obtained as a ^ 0, is described by the 
equation 

Lue = in D, L :=V^ + k^n^{x) ■.= V^ + - q{x), (22) 
n'^{x) := nQ{x) — k'''^p{x), p{x) := 47rA^(x)/i(x), q{x) := qo{x) +^(2;), (23) 

Ue = u{x) + Veix), (24) 

where Ve{x) := jj^g{x,y)p{y)ue{y)dy. 

The scattering amplitude, corresponding to n^{x), (or, which is equivalent to the po- 
tential q{x) ) is: 



A{p, a) = Ao(/?, a) + ^ u{y, -p)p{y)ueiy)dy, (25) 



47r 



where u{x, —[3) is the solution to problem ((T])-([3]) with uq = , that is, the scattering 

solution in the absence of the embedded particles when the incident direction of the plane 
wave Uq equals to —(3. 
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Remark 1 Formula (|25p is based on the Ramm's lemma (see [21 p. 46]) which says 



^ik\y\ I y 

9{x, y) = —^Mx, -P) + 0(--2 ), \y\ ^ oo, /3 := — , (26) 
4vr|y| |y|^ \y\ 



where u{x, —(5) is the scattering solution with the incident direction —[3 € S"^ and the 
remainder O(j^) is uniform with respect to x running through any fixed bounded domain. 

We can now formulate a recipe for creating material with a desired refraction coefficient: 
Recipe: To obtain material with a desired refraction coefficient n^(x), given 
a bounded domain D, filled with a material with the coefficient no(x), one 
embeds into D small balls Dmi centered at Xm and of radius a, distributed 
according to ([5]), with 4>{a) defined in (fT2]l . so that conditions (i9l)- (fTT]) hold. 
Then one finds h{x) and N[x) from the equation 

k^[nl{x) - n^{x)] := p{x) = iTTh{x)N{x). (27) 

Equation (|27p for N{x) > and h{x),lmh < 0, where p{x) is given, has infinitely 
many solutions. For example, one can choose N(x) > arbitrary, and then find 
uniquely hi and /12 by the formulas: 

hi ■.= Reh = ^ ^l. , , /i2 := Im/i = — (28) 
47riV(x) 47riV(x) ^ ^ 

where p = Pi + ip2i Pi = Rep, P2 = Imp. 

By Theorem [5] the material with the embedded small particles has the 
desired coefficient n^(x) with the error that goes to zero as a —5^ 0. 

We do not give proofs of Theorem [THSj These proofs can be found in the cited papers. 
Our goal is to formulate clearly the recipe for creating materials with a desired n^(x), so 
that engineers and physicists can try to implement this recipe practically. 



3 A discussion of the recipe 

There are two technological problems that should be solved in order that the recipe be 
implemented practically. 

Problem 1: How does one embed many small particles in a given material so that the 
desired distribution law (j5|). (|lip - (|12p is satisfied? 

Problem 2: How does one prepare a small particle with the desired boundary impedance 

/- h{Xm.) ip 

(,m — • 

The first problem, possibly, can be solved by stereolitography. The second problem 
one should be able to solve because the limiting cases Cm = (hard particles) and Cm = 00 
(soft particles) can be prepared, so that any intermediate value of Cm one should be able 
to prepare as well. The author formulates the above technological problems in the hope 
that engineers get interested and solve them practically. 
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3.1 Negative refraction 

Material with negative refraction is, by definition, a material in which group velocity is 
directed opposite to the phase velocity, (see [Ij and references therein). Group velocity 
is defined by the formula Vg = Vkw(k). Phase velocity Vp is directed along the wave 
vector k° = In an isotropic material lo = w(|k|), and lo = ;^^^|jy, so u;n(x,u;) = c|k|. 
Differentiating this equation yields 

VkW [n{x, uo) + = ck°. (29) 

OLU 

Thus, 

Vflh + o;— ] =ckO. (30) 

Wave speed in the material with refraction coefficient n{x,uj) is |vp| = where c is 

the wave speed in vacuum (in D'). 

For Vg to be directed opposite to k*^, that is, opposite to Vp, it is necessary and sufficient 
that 

c)ti 

n + uj—<0. (31) 

OUJ 

If the new material has n{x,uj) satisfying (j3ip . then the new material has negative refrac- 
tion. 

One can create material with n{x,uj) satisfying pip by choosing h = h{x,uj) properly. 
Namely, (I27p implies 

n^{x, io) = nl{x) - A7:k-^Nix)h{x, lo). (32) 
Assuming Uq > A7rk^'^N{x)h{x,io), Imh = 0, one has ()3ip satisfied if 

/ 27rA;~2iVfxl^^^^ 

^nl-A7rk-mix)hix,u;) < , ^ \' />^ ==, 

V ~ 4:TTk-^N{x)h{x, uj) 

or 

, , n?.(x) dh , , , 

This inequality has many solutions, e.g., the function h{x,uj) = e^^'^/io + ^e^'^t(x) solves 

Let us formulate the technological problem solving of which allows one to implement 
practically our method for creating materials with negative refraction. 

Problem 3: How does one prepare a small particle with the impedance Cm = '^^^^'^^ , 
where h{x,Lo) satisfies ([33]) ? 
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